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Q ABSTRACT 

Context. We present an analysis of the kinematics and physical properties of the Class I driven jet 
HH46-47 based on IR medium and low resolution spectroscopy obtained with ISAAC on VLT. 
Aims. Our aim is to study the gas physics as a function of the velocity and distance from the 
source and compare the results with similar studies performed on other Class I and classical T 



> 
CO 

Tauri jets, as well as with existing models for the jet formation and excitation. 
Methods. The ratios and luminosities of several important diagnostic lines (e.g. [Fen] 1.644, 



CN 
CN 

o 



1.600pm, [Pn] 1.189yum, and H 2 lines) have been used to derive physical parameters such as 
electron density, H 2 temperature, iron gas-phase abundance and mass flux. [Fen] 1.644yum and 
H 2 2. 122/im Position Velocity Diagrams (PVDs) have been in addition constructed to study the 
kinematics of both the atomic and molecular gas. 

Results. Within 1000-2000 AU from the source the atomic gas presents a wide range of radial 
velocities, from ~-230kms~' to ~ lOOkms -1 . Only the gas component at the highest velocity 
(High Velocity Component, HVC) survives at large distances. The H 2 , shows, close to the source, 
only a single velocity component at almost zero velocity, while it reaches higer velocities (up to 
~95 km s ) further downstream. Electron densities (n e ) and mass ejection fluxes (M je( ) have been 
measured separately for the HVC and for the component at lower velocity (LVC) from the [Fe n] 
lines. n f increases with decreasing velocities with an average value of ~ 6000 cirr 3 for the LVC 
and ~ 4000 cirr 3 for the HVC, while the opposite occurs for M ie , which is ~ 0.5-2x 10~ 7 M yr _1 
and ~ 0.5-3. 6x 10~ 8 M yr~' for the HVC and LVC, respectively. The mass flux carried out by the 
molecular component, measured from the H 2 lines flux, is ~ 4 x 10~ 9 M yr~' . We have estimated 
that the Fe gas phase abundance is significantly lower than the solar value, with ~ 88% of iron 
still depleted onto dust grains in the internal jet region. This fraction decreases to ~58%, in the 
external knots. 

Conclusions. Many of the derived properties of the HH46-47 jet are common to jets from YSOs 
in different evolutionary states. The derived densities and mass flux values are typical of Class 
I objects or very active T Tauri stars. However, the spatial extent of the LVC and the velocity 
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dependence of the electron density have been so far observed only in another Class I jet, the 
HH34 jet, and are not explained by the current models of jet launching. 

Key words, stars: circumstellar matter - Infrared: ISM - ISM: Herbig-Haro objects - ISM: jets 
and outflows - ISM:individual objects: HH46-47 



1. Introduction 



Protostellar jets are tightly associated with the first stage of protostar evolution. They are believed to 
remove angular momentum, disrupt infalling material from the cloud, inject turbulence in the ISM 
and modify its chemistry. Little is known, however, on how protostellar jets originate. There are two 
principal schools of tho ught regarding where protostel lar jets are launched from: the X-wind and 



the Disc-wind models dShu et al 



1995 



Ferreira 



1997). While disc-wind models assume ejection 



from a large spread of radii from the disc, X-wind models assume that the ejection of material takes 
place only at a single annulus at the interaction region between the stellar magnetosphere and the 
inner edge of the disc. Both models were born in an attempt to explain the kinematic al and physical 
properties observed at op tical wavelengths in jets from classical TTauri stars (e.g. 



Bacciotti et al 



2000; 



Lavallev et al. 



1997ft . Very few studies have been done, however, regarding the physical 
properties of Class I jets. Since these objects are still highly embedded in their molecular clouds, 
they are affected by a large extinction. Thus, it is not possible to detect the region near to the exciting 
source at optical wavelengths. In this context, infrared spectroscopy is a very useful tool to probe 
the physical properties and kinematics of Class I jets close to their base, where their characteristics 
have not yet been modified by the interaction with the surrounding material. 

Several studies have beei 
tic on protostellar jets (e.g., 



Nisini et al. 


2005; 


Podio et al. 


2006; 


Davis et al. 


2006: 


Takami et al. 



2006). The main lines used to probe physical and kinematical properties in NIR are [Fen] (e.g. 
1.644 /mi, 1.600 pun) and H2 lines (e.g. 2. 122 /mi). These studies revealed that the kinematical 
structure at the jet base of Class I and classical T Tauri stars (CTTSs) jets is very similar. Both 
show a broad line emission near the source consisting of a collimated and large scale jet at high 
velocity (the so-called HVC) and a slower velocity component (the so-called LVC) that is only 
detected around the central source. Noticeably, this kind of velocity structure is often observed 
in both the atomic and molecular components of the jets, th e latter probed thro ugh the H2 emis- 
sion and called Molecular Hydrogen Emission Line regions dDavis et al.ll2001bl) . Concerning the 
physical properties, Class I jets show higher densities and mass loss rates than CTTS jets, as ex- 
pected for less evolved objects. Both disc-wind and X-wind models predict the presence of a large 
spread of velocity at the j et base as shown by few synthetic pos ition velocity diagrams (PVDs) con- 



structed for both models dShang et al 



1998 



Pesenti et al 



2004), although attempts to reproduce the 



specific characteristics (e.g. kinematics, presence of the MHEL region, physical parameters) dis- 
played by the Class I jets are still lacking . We have now started a project aimed to determine few 
basic physical parameters of embedded jets as a function of both velocity and distance from the 
driving source, using near-IR (NIR) medium resolution (R~10000) spectra taken with ISAAC on 
the Very Large Telescope (VLT). The goals of this study are to make a quantitative comparison 
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* Based on observations collected at the European Southern Observatory, Paranal, Chile (ESO programmes 
0.74.C-0286(A)). 
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Fig. 1. Position of the ISAAC slit superimposed on the H2 2.12/mi image of the HH46-47 jet. At 
the source (HH46-IRS) and knot S position a negative residual from the continuum substraction is 
present. In dividual knots along the jet are identified following and our own (see text) nomenclature 
and that of Eisl offel & Mundt dl994l) . The northeastern knots correspond to the blue-shifted jet. 



between the physical properties of Class I and CTTS jets and to constrain the origin of the diffe- 
rent velocity components and the mechanism responsible for their excitation. In a previous paper 
( Garcia Lopez et al]|2008 ). we have applied our NIR diagnostic techniques to two jets, HH34 and 
HH1. 

We now extend this previous study to the HH46-47 jet. This known Class I jet is at a distance 
of ~450pc emerging from a dense reflecting nebula. The jet is located at the border of a Bok 
globule and powers an atom ic and molecular bip olar flow. The driving source (HH46-IRS) is in 
a phase of high accreti on ( Antoniucci et al.ll2008l) and is known to be a deeply embedded binary 
system dReipurth et alj 12000). The northeastern (blue- shifted) portion of the flow is detected at 



optical wavelengths giving rise to strong HH objects (IHartigan et al 



1993 



2005 ). Recently, the 



Nishikawa et al 



kinem atics of the internal region of this optical jet has been studied in details by 
(2008). At infrared wavelengths, the HH 46 jet and its red-shifted counter-jet are detected down to 
the central source (lEisloffel et al.l ll994). thus this object represents a very suitable target for our 
NIR spectroscopy study. 

In Sect. 2, we present the observations and the data reduction. In Sect. 3 we describe the results 
of our observations, while the kinematics of the atomic and molecular components are presented 
in Sect. 4. In Sect. 5 and 6 we describe the spectral diagnostic analysis applied to ISAAC spec- 
tra obtained at low and medium resolution, respectively. We have compared the results found on 
the HH46-47 jet with the HH34 jet in Sect. 7 and the main conclusions of our work are finally 
summarised in Sect. 8. 
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Table 1. Lines observed in the knot Zl of the HH46-47 jet. 



Line id. 


A 


F 


AF 


Line id. 


A 


F 


AF 




(jum) 


(io- 16 


erg cirT 2 s~') 




(jum) 


(io- 16 


erg crrT 2 


[Till] 2 G7/2" 4 F7/2 


1.140 


8.3 


3.5 


H 2 1-0 S(6) 


1.788 


9.4 


0.5 


H 2 2-0 S(l) 


1.162 


5.8 


2.2 


[Fen] a 4 D 3 / 2 -a 4 F 3 / 2 


1.797 


9.8 


1.1 


[PuJ^-'Dz 


1.188 


3.5 


1.0 


[Fen] a 4 D 5 / 2 -a 4 F 5 /2 


1.800 


7.1 


0.8 


[Fen]a 4 D 7/ 2-a 6 D 9/ 2 


1.257 


106.0 


0.8 


[Fen] a 4 D 7/2 -a 4 F7/2 


1.810 


38.2 


2.1 


[Feii]a 4 D 1/ 2-a 6 D 1/2 


1.271 


5.9 


1.0 


H 2 2-1 S(5) 


1.945 


11.0 


1.5 


[Feii]a 4 D 3/ 2-a 6 D 3/ 2 


1.279 


7.8 


0.1 


H 2 2-1 S(4) 


2.004 


2.7 


0.3 


[Fen]a 4 D5 / 2-a 6 D 5/ 2 


1.295 


11.9 


0.9 


H 2 1-0 S(2) 


2.034 


26.1 


0.6 


[Feii]a 4 D 3/ 2-a 6 D 1/2 


1.298 


6.7 


0.1 


H 2 2-1 S(3) 


2.074 


8.1 


0.5 


[Fen] a 4 D7/2-a 6 D7/2 

L J / 1 Z, / 1 Z, 


1.321 


306.0 


0.8 


H 2 1-0 S(l) 


2.122 


68.5 


0.3 


[Fenl a 4 D5/2-a 6 D 3 /2 


1.328 


8.4 


0.9 


H 2 2-1 S(2) 


2.154 


3.9 


0.3 


[Fen] a 4 D 7 /2-a 6 D5/ 2 


1.372 


71.7 


1.0 


H 2 2-1 S(l) 


2.248 


8.5 


0.5 


[Fen] a 4 D5/2-a 4 Fa/2 


1.534 


173.0 


0.7 


H 2 3-2 S(l) 


2.386 


1.6 


0.4 


[Fen] a 4 D3/2-a 4 F 7 /2 


1.600 


10.6 


0.5 


H 2 1-0 Q(l) 


2.407 


50.0 


0.8 


[Fen] a 4 D 7/ 2-a 4 F 9 /2 


1.644 


160.0 


0.6 


H 2 1-0 Q(2) 


2.413 


18.3 


0.9 


[Fen] a 4 Di /2 -a 4 F 5/2 


1.664 


7.0 


0.7 


H 2 1-0 Q(3) 


2.424 


47.4 


0.9 


[Fen] a 4 D 5/ 2-a 4 F 7/ 2 


1.677 


142.0 


0.3 


H 2 1-0 Q(4) 


2.437 


18.9 


1.3 


H 2 1-0 S(9) 


1.687 


4.2 


0.5 


H 2 1-0 Q(5) 


2.455 


42.9 


1.5 


[Fe n] a 4 D 3/ 2-a 4 F 5/2 


1.712 


2.1 


0.4 


H 2 1-0 Q(5) 


2.455 


42.9 


1.5 


[Fen]a 4 Di /2 -a 4 F* /2 


1.745 


6.3 


0.8 


H 2 1-0 Q(7) 


2.500 


47.4 


3.7 


H 2 1-0 S(7)* 


1.748 


150.0 


0.6 











*Blend of the two indicated lines. 
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Fig. 2. LR spectrum from 1.00 to 2.50 /urn of the knot Zl in the HH46-47 jet. The stronger lines 
are identified. 
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The observations were made on 29-30 December 2004 at ESO VLT telescope on Cerro Paranal, 
Chile, using the infrared spectrograph and camera ISAAC. The spatial scale of the instrument is 
07 148/pixel. We have taken low resolution spectra in the J-band (1. 1-1.4 fxm) and both, medium 
and low resolution spectra in the H ( 1.42-1. 82 yum) and K (1.82-2. 5/im) bands. The adopted slit 
width were 076 and 073 for the low and medium resolution, respectively, corresponding to nominal 
resolutions of 860 (in J), 840 (in H) and 750 (in K) at low resolution (LR) and 10000 (in H) and 
8900 (in K) at medium resolution (MR). The slits were all aligned along the jet, with a position 
angle (RA.) of ~57°. Total integration times were 600 s for each filter in the LR mode and, 3600 s 
and 7200 s for the K- and H-band in the MR mode. 

Data reduction was performed using standard IRAb3 tasks. Wavelength calibration was done 
using the atmospheric OH emission at MR and an argon lamp at LR. The atmospheric spectral 
response was corrected by dividing the object spectra by the spectrum of a telluric standard star for 
both modes. Flux calibrations were performed using photometric standards observed at an air mass 
and seeing similar to our target observations. This procedure probably leads to some flux losses due 
to the finite slit width with respect to the jet width, especially in the case of the bow shaped knot 
Zl. In MR spectroscopy, the continuum of the HH46-IRS source was subtracted using the IRAF 
task BACKGROUND. 

The reduced MR spectra on the source position have been already presented in 



Antoniucci et al 



(2008) (reduced and analysed in the same way as explained above) and thus they will not be dis- 



cussed in the present paper. 



3. Results 

Fig. Q] shows a H2 2.12/mi continuum-subtracted image of the HH46-47 jet. The slit position 
is depicted over the ima ge and the covered knots are named following the nomenclature of 



Eisloffel & Mundt ( 1994). In addition, we name as knot S the emission that peaks at the source 
position and extends for ~3" in the blue-shifted jet. The observations cover the blue-shifted jet 
knots HH46-1 and 2 and the knot B8, which is l ocated at the apex of a bo w shock. In the red- 



Eisloffel & Mundt ( 1994) that have been here 



shifted lobe, the slit comprises knots Z and Y in 
separated in knots Zl and Z2, and Yl, Y2 and Y3. 

The low resolution spectra present several lines of both atomic and molecular species in the 
different knots. Figure|2] shows the complete spectrum of knot Zl, which is among the brightest 
ones along the jet. A list of the observed lines, fluxes and identification in knot Zl is also reported 
in TableQ] while Tables |A~T| and |A.2"l in the AppendixlAllist all the lines detected in the other knots. 
[Fen] transitions from the first 12 fine-structure levels, and H2 transitions from the two first vibra- 
tional levels have been detected in all of the knots at various intensity levels. The knot S is, together 
with the knot Zl , the richest in ionic lines. Here [P 11] 3 P2 - 1 D2 1.1 89 yumhas also been detected. 

The medium resolution K-band spectra show several lines, the brightest being H2 1-0S(1), 3- 
2S(4) and 2-lS(2), [Fe 11] 2. 1349 and 2. 1609 yum and Bry (see, Fig.0. In addition, the [Fe n] 1 .600 



1 IRAF (Image Reduction and Analysis Facility) is distributed by the National Optical Astronomy 
Observatories, which are operated by AURA, Inc., cooperative agreement with the National Science 
Foundation. 



._. 
o 

CD 
W 

o 

CD 
O 

|-10 

V) 
CD 

I" 15 

o 

g "20 
§■25 



Garcia Lopez, R. et al.: Velocity resolved observations of HH46-47 

HH 46-47 jet, K-band 
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Fig. 3. MR continuum subtracted spectral image of the HH46-47 jet in the K-band. The distance 
from the source is represented on the left-hand Y-axis, while the knots are labelled on the right- 
hand. 



21600 



21700 



and 1.644 /an lines have been detected in the H-band spectra. The presence of the Bry near the 
source and the H2 3-2S(4) line at a distance of 10" to 25" from HH46-IRS testifies to the high 
excitation conditions in these regions. In additio n, the on-source spectru m shows other features, as 
the Brl4 and Brl2 and CO transitions (see also, Antoniucci et al.ll2008 ). 



4. Kinematics 

In order to study the kinematics of the atomic and molecular component of the HH46-47 jet as 
a function of distance from the central source, we have constructed Position Velocity Diagrams 
(PVDs) of the [Fen] 1.644yum and H2 2.122/im lines (see, Fig.|5j from the MR spectral images. 
The velocity has been expressed with respe ct to the Local Stan dard of Rest (LSR) and corrected 



for a parental cloud velocity of +20 km s 1 (IHartigan et al 



1993b . In the Y-axis the distance from 



the driving source HH46-IRS, in arcsec, is represented. The continuum of the source has been 
subtracted in both the PVDs to more clearly see the structure of the lines near HH46-IRS. The 
intensity scale is different in the two PVDs in order to show the different condensations in the 
atomic and molecular gas. 

Emission down to the central source is detected for both the atomic and molecular components 
(knot S). The [Fen] emission broaden as it approaches the source position and spans velocities 
from ~ -200 km s _1 down to km s _1 within 2 " from the source. Inside 1 " there is also emission at 
red-shifted velocities up to ~ +100kms~'. We will call the blue-shifted component at high velocity 
that can be traced up to large scale the high velocity component (HVC) and the more compact 
emission component from ~-150kms 1 to ~0kms _1 the low velocit y component (LVC). These 



two components have been also evidenced in Her optical observations (INishikawa et al 



2008). At 
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Knot 


[Fen] 1.644 /an 


H2 2.122/im 




v(kms ) 


v(kms _1 ) 


B8 


-177 




HH46-1 


-235 


-20 


S 


-213,(-138) 


-15 


Z2 


+ 144, +92 


+26 


Zl 


+ 112 


+89, +10 


X 


+94 


+93, +14 


Y3 




+94 


Y2 




+92 


Yl 




+92 



Note: Radial velocities corrected from a cloud velocity of +20 km s 1 (lHartigan et al.llT993h and 
with respect to the LSR. The radial velocity error is 2 km s _1 . 



variance with the ionic emission, the H2 emission close to the source (knot S) shows only one 
component at nearly km s _1 . 

The radial velocities of the [Fen] 1.644/im and H2 2.122/vm lines are reported in Table|2] They 
have been computed by a Gaussian-fit to the line profile of every knot. The velocities of the [Fe n] 
emission range from -235kms~' to -177 kms~' and from +144 kms -1 to +94kms _1 in the blue- 
and red-lobes, respectively: thus in both lobes the peak velocity decreases as the distance from the 
driving source increases in both lobes. The HVC and LVC in knot S peak at -213 km s _1 and around 
-138 kms -1 . At knot Z2 the emission splits in two separate velocity components (at +144 kms" 1 
and +92kms~ 1 ) whose origin is likely connected with the presence of the large bow shock struc- 
ture. 

The radial velocities of the H2 emission vary from -15 km s _I to -48 km s in the blue-lobe an 
from +26kms~ 1 to +94kms~ 1 in the red-lobe. At variance with the [Fen] emission, knot S only 
shows a single velocity component at -15kms _1 while both knots Zl and X present two velocity 
components at ~ +90kms~' and +10kms~'. The high and low velocity components present in 
these knots and the broad profiles seen in H 2 (see, Fig.[4]l are similar to the profiles observed in 



other molecular jets (e.g. jDavis et alj|2001a[) . As shown in Fig.|4l the relative contribution of the 



high and low velocities components to the total flux changes from knot Zl (where the component 
at Okms -1 is larger than that at high velocity) to knot X (where the component at 0kms~' is 
significantly reduced). These profile s are similar to tho se expected from b ow-shocks structures 
viewed at a certain inclination angle ( Schultz et al.ll2005 ). With respect to the Schultz et al. (2005) 



models, the HVC and LVC contribution to the total flux depends, in our case, not only on the 
viewing angle of the bow-shock but also on how the H2 cooling is distributed from the apex to the 
wings of the bow. In particular, the profiles of the knots Zl and X are consistent with bow-shocks 
seen at 30° with respect to the line of sight, but with a different relative contribution to the H2 
cooling between the apex and the wing. In knot X the apex mostly contribute to the total flux while 
in knot Zl the situation is the opposite. Such an effect could be actually due to a not good alignment 
of the slit with the bow-shock apex in knot Zl. 
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Fig. 4. Normalised line profiles of the H2 2. 122 /im line for knots Zl (upper panel) and X (lower 
panel). In the X-axis the LS R radial velocity is r epresented. The velocities are corrected from a 
cloud velocity of +20kms _1 (Hartiga n et alJI 1993b 



Table 3. A v along the HH46-47 jet. 



Knot 


Ay ± AAi£ 




(mag) 


S 


6.6+0.2 


Z2 


12.5+1.6 


Zl 


6.1+0.1 



" Visual extinction measured from the ratio [Fen] 1.64/1.25 jim 



5. Low resolution spectral analysis 

5.1. Extinction 

Visual extinction values can be derived from the [Fen] 1.644/1. 251 ratio: these lines share the 
same upper level and thus, their theoretical ratio is a function of the frequencies and probabilities 
of the transition and are independent of the plasma conditions. To conve rt the observed ratio s into 
A v values, we have adopted the radiative transition probabilities given bv lOuinet et al.l ([1996) using 
the relativistic Hartree-Fock (HFR) computer program. The extinction errors have been estimated 
directly from the flux error of each line. The results are reported in Table[3]for the individual knots 

where the S/N of the two lines is larger than 10. 

Ay varies from ~6 in the knots Zl and S to ~ 12 in the knot Z2. iFernandesI (120001) found a 



value of Ay = 9.38 + 1.49mag in a region that covers our knots HH46-1 and S. This value has 
been ho wever derived employing the theoretical [Fen] 1.644/1 .257 /an ratio given by th e radiative 



rates of 



Nussbaumer & Storevl 019881) that predict higher values with respect to the 



Ouinet et al 
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Fig. 5. Continuum-subtracted PVDs for the [Fen] 1.644 fan and H2 2.122 /an lines along the 
HH46-47jet. Contours show values of 3, 9, . . . , 243 a for the [Fen] line and 4.5, 13.5,. .. ,364.5 a 
for the H2 . On the Y axis distance from the source HH46-IRS in arcsec is reported. 



( 1996) rates (iNisini e t al. 2005). On the other hand 



Antoniucci et al 



(2008) found an A v -40 mag 



on-source, indicating a large density gradient towar ds the central object. One should note here 
that the extinction found by 



Antoniucci et al 



(2008) on the basis of photometric measurements 
coupled with Bry emission, is larger than the one we found on the knot S, since [Fen] lines trace 
more external regions. In addition, the knot S covers a region that extends ~3" around the source 
position. 
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Table 4. dv„ T, N(H 2 ), M Hl along the HH46-47 jet. 



Knot 



Zl 
X 
Y3 
Y2 
Yl 



" Tangential velocities from 
ejection rate 



(kms- 1 ) 



jb 

(K) 



N(H 2 ) C M H / 
(10 17 ctrT 2 ) (Moyr" 1 ) 



131 
19 
122 
122 
122 



2797 + 100 
2078 ± 42 
2564 + 53 
2330 + 50 
2093 + 42 



3.3 
3.7 
3.6 
5.1 
2.4 



4.5 xlO" 9 
7.3xl0~ 10 
4.6xl0~ 9 
6.4x1 0- 9 
3.0xl0~ 9 



Eisloffel & Mundt ( 1994). H 2 temperature. r H 2 column density. H 2 mass 



The extinction in the knot Z2 is roughly twice the extinction of the knots S and Zl. This 
fact could be due to the collision of the jet with stationary material on a side of the cavity. 

( 1996) note the presence of an oblique shock due to the impact of the blue-shifted 



Heathcote et al 



lobe of the jet with the cavity at ~6" NE of the source. The same process may occur at the same 
distance in the red-shifted counter-lobe, giving rise to an increase in extinction. 



5.2. H 2 temperature and mass flux 

Using the different H2 detected transitions, excitation diagrams for several knots (Zl, X, Y3, Y2 
and Yl) along the jet showing enough H2 lines have been constructed. As an example, Fig.|6]shows 
the excitation diagram of the knot Zl, the brightest knot observed in HH46-47. The excitation dia- 
grams for the rest of the knots are reported in the AppendixlBl In this kind of diagrams the natural 
logarithm ofN(v,J)/ g v j versus the excitation energy E( v , J) of a given level (v , J) is represented, 
where N(v , J)/g v j is the ratio of the column density over the statistical weight. The different sym- 
bols indicate lines coming from different vibrational levels, while the straight line represents the 
best linear fit through our data. The slope of the line is proportional to the inverse of the gas tem- 
perature. Extinction values as derived from the [Fe 11] lines have been adopted. The extinction can, 
in principle, be also derived from the H2 1-0S(1)/1-0Q(3) line ratio. However, this method involves 
lines close in wavelength and thus it is not sensitive enough to extinction variations. In addition, the 
Q lines fall in a poor atmospheric transition region (around ~2.3 /mi) being their fluxes affected by 
large uncertainties. The fact that our data can be well traced by a straight line on the ro vibrational 
diagram indicates, however, that the extinction values derived from the [Fe n] lines are in good 
agreement with the H2 ones. The resulting temperature values are indicated in Table 4 and range 
between -2000 K and -2800 K. 

The mass flux carried out by the warm molecular component (M//,) has been also estimated 
from the observed H2 transitions in the knots Zl, X, Y3, Y2 and Yl. can be inferred from 
the relationship M# 2 =2 N(H2) A dv t /dl,. Here, A is the area of the emitting region sampled by the 
slit, dl, and dv, are the projected length and the tangential velocity, and Nify) is the total column 
density. M#, has been derive d in the same spatial regi on extracted for the [Fe 11] line analysis. The 
dv, value has been taken from Eisloffel & Mundt ( 1994), while N(H2) for each knot has been found 



from the intercept to zero of the straight line fitted to the H2 transitions in the Boltzmann diagram. 
The derived values for each knot are reported in Table[4] The average M#, value along the jet is 
~ 3.8 x lO _9 M yr _1 . This estimate is however a lower limit since we are only considering the 
warm component of the gas. 
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Knot Zl 
T = 2797±100 K A„ = 6 maa 




E V>J (10 3 K) 

Fig. 6. H2 rotational diagram for the different lines in the knot Zl of HH46-47. Different sym- 
bols indicate lines coming from different vibrational levels: circles indicate v=l, triangles v=2 and 
squares v=3. The straight line represents the best fit through the data. The corresponding tempera- 
ture and the adopted extinction value is indicated. 



Table 5. Electron density along the HH46-47 jet. 



Knot 


(") 


n e ± An^ 

(cm" 3 ) 


n e (HVC)^ 
(cm 4 )" 


n e (LVC)* 
(cm" 3 ) 


S 


(-3.4,-1-2.8) 


6300 +170 ° 

"•'""-1300 


5400 


7650 


Z2 


(-9.3,-5) 


1300 +180 ° 

1J ""-1200 


2350 


4600 


Zl 


(-14.8,-9.3) 


3200 +70 ° 

JZ - UU _850 


3100 





" Distance from the source in arcsec. Negative values correspond to the southeastern, redshifted jet axis. 

* Electron density measured from the ratio [Fen] 1.60/1.64 yum . 

Electron density measured from the ratio [Fen] 1.60/1. 64 /an for the HVC and LVC. 



5.3. Electron density 



The electron density (n c ) can be derived from [Fen] lines coming from levels with similar exci- 
tation energies. This is the case of the ratio [Fen] 1. 600/1. 644/vm, which is sensitive to n e values 
between 10 3 - 10 5 cra- 3 and weakly dependent on the temperature . This ratio has been fitted with 
a NLTE code including the first 16 fine-structure levels of the [Fen] ( Nisini et alj|2002 ). 

The electron density has been derived integrating the total flux of the line profile of every 
extracted knot for both lines. The n e values and the intervals used to extract the knots are reported 
in Table|5] Previous studies on protostellar jets, likeHHl, H H1 1 1 and HH34, have shown evidences 
of a decrease of the elect ron density with distance (see, e.g 



Bacciotti & Eisloffel 1999; 



Nisini et al 



2005 



Podio etal 



2006). The same behaviour may also occur in HH46-47. Indeed the electron 
density roughly decreases from a value of 6300cm- 3 in knot S to a value of 3200cm 3 in knot Zl. 
These two knots are displaced, however, in different lobes of the jet. In addition, the increase of 
n e from knot Z2 to knot Zl and the the lack of further [Fe 11] emitting knots prevent us to properly 
probe the dependence on distance of n e . 
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Knot S 




0.0 

-300 



-200 -100 
velocity(kn / s) 

Knot Z2 




100 200 
velocity(kn / s) 

Knot Z1 



100 200 
velocity(km / s) 



.500 




300 



Fig. 7. Normalised to one line profiles (lower panels) of the [Fen] lines 1.644yum (solid line) 
and 1.600/im (dotted line), and the electron density in each velocity channel (upper panels) for 
different extracted knots along the HH46-47 jet. 



5.4. [Fen] depletion 



Refractory such species as iron are expected to be locked into dust grains in the quiescent conditions 
of molecular clouds. High velocity shocks inside the jet beam are, however, supposed to destroy 
the grains and release Fe atoms in gaseous form. In recent years several studies have been done in 
order to quantify the presence of dust along protostellar jets and measure the shock efficiency i n 
disrupting the dust grain cores (see e.g. 



Nisini et al. 2005 



Podio et al. 2006; 



Giannini et al. 2008). 



The iron gas phase abundance can be measured comparing the emission of [Fe n] with that of a non- 
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refractory species of known abund ance. A very su itable ratio is the one between the [Fe n] 1 .257 /mi 
line and the [Pn] 1.1 88 /mi line dOliva et al.ll2.Q0lh . These two lines are very close in wavelength 
and thus their ratio is poorly affected by extinction. In addition, both lines are excited under similar 
conditions and can be assumed in the first ionised state at similar percentage. Assuming a solar 
abundance for bo th elements, the expected theoretical line ratio is [Fen]/[Pn]~l/2 [Fe/P]~56 



dOliva et al 



120011) . Thus, observed values lower than this indicate that some of the iron is locked 
in grains. 

Along the HH46-47 jet, both the [Fen] 1.257 /an and the [Pn] 1.1 88 /mi lines have been de- 
tected in knots S and Zl. The measured [Fen] / [Pn] ratio is ~7.3 in knot S and ~25.6 in knot Zl. 
These values indicate that around 88% and 58% of iron is still in dust grains in knots S and Zl, re- 
spectively. This leads to a ratio (Fe/H) gas ~ 3.92 x 10~ 6 in the knot S and a (Fe/H) gflJ ~ 1.18 x 10~ 5 
in the knot Zl. The str ong iron depletion at the jet base is still lower than the 99% diffuse gas 
iron depletion (see, e.g., Savage & Sembachll996 ). indicating that a small fraction of the iron have 
come back into gas phase at the base of the jet with respect to the interstellar value. Evidence that 
the gas phase iron abundance is low in knots close to the source and increases with the distance 
have been also found in the HH1 jet ( Nisini et alj|2005 ). 

In knot Zl, the [Tin] 1.140//m line has been also detected. Ti is expected to be fully ionised 
in the jet plasma, since it has a low ionisation potential of 6.82 eV. As iron, the titanium is a high- 
refactory element and thus the Fe/Ti ratio can give us some clues about the relative abundance of 
both elements, as well as to the selective depletion pattern of both elements in the dust along pro- 
tostellar jets. With this aim we have comp uted the [Fe n] 1 .257 /an / [Tin] 1.140//m intensity ratio 
using a NLTE model for these species dGarcia Lopez et al.l 120081) . The observed ratio, assuming 
the physical conditions present in knot Zl is ~70, indicating a Fe/Ti ra tio 5 times smaller th an the 
theorically expected of [Fe/Ti] = 354 using the solar abundances of Asplundet al. (2005). The 
small Fe/Ti is consistent with what found in the HH1 jet ( Garcia Lopez et all2008 ). where a Fe/Ti 
ratio 2-3 times lower than solar was derived. 



6. Medium resolution spectral analysis 

6. 1 . Electron density and Mass flux 

The [Fen] 1.600 and 1.644 /mi lines have been also detected at medium resolution in the knots 
S, Zl and Z2. In these knots, we can test how the electron density varies as a function of the 
radial velocity, by measuring the 1.600/1. 644 /mi ratio in each pixel along the spectral profile. 
The [Fe n] ratio was computed only for the velocity points where the intensity is greater than 
S/N ~ 3. Figure|7] shows electron density variations along the MR spectral profile of the knots 
S, Zl and Z2 (upper panels) together with the normalised line profile of the [Fen] 1.600//m and 
1.644//m lines (lower panels). The electron density increases as the velocity decreases in all of the 
knots where two velocity components were detected in [Fe n]. In particular, the knot S close to the 



source has a behaviour very similar to what was found in HH34 dGarcia Lopez et al 



2008), with 



the electron densities higher in the gas at low velocity. We have separately measured the n e in the 
HVC and LVC components. Near to the source (knot S), we have defined a HVC ranging from 
~-239kms~' to —181 kms~' (consistent with the FWHM of a Gaussian-fit through the line), and 

-50 km s _1 . Far from the source, only the knot Z2 shows a clear 
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Table 6. M jet along the HH46-47 jet. 



Knot 


dv t (HVC^ 




dv r (LVC) fl 






km s _I 


1(T 7 M yr 1 


kms 1 


lO^Moyi- 1 


S 


316 


0.3-1.9 


205 


0.5-3.6 


Z2 


213 


0.4-0.9 


136 


0.2-0.4 


Zl 


166 


0.4-1.0 







" Tangential velocity derived from the observed radial velocities assuming an orientation angle of the jet 
with respect to the plane of the sky of 34° (Eisloffel & Mundt 1994|). 

* Mj e , measured from the luminosity of the [Fen] 1.64/jm for the HVC and LVC. The lower value has 
been derived assuming an iron solar abundance, while the upper value has been inferred assuming the iron 
gas-phase abundance derived from the [Fe n]/[P n] ratio. 



double velocity component: here we have defined the intervals from ~110kms-' to ~167kms _1 
and from ~80km s _1 to ~1 10 km s _1 for the high and low velocity gas, respectively. 

The n e values for the integrated velocity components are reported in Table|5] The elec- 
tron density of the LVC (n c (LVC)~4600-7600) is larger than the electron density of the HVC 
(n c (HVC)~2300-5000) as evidenced by Fig. [7] The behaviour found in knot S is similar to what we 
have already found at the HH34 jet base, and contrast with the results on T Tauri stars, where the 
highest electron density is associated with the high velocity gas pertaining to the more collimated 
inner jet section, w hile the low velocity component is a less collimated gas at lower density (e.g., 



CofFev et al. 



2008). These contrasting results may reflect a different excitation conditions in the 



structure of Class I jets, or may be a result biased by the low spatial resolution of our observations 
that does not allow us to resolve the same spatial scales probed in T Tauri jets, namely within 
1 00-200 AU. 

At the bow shock positions far from the source, like in the Z2 knot, variations of the electron 
density with velocity ma y depends on the bow inclination angle and section of the bow intercepted 
by the slit, as shown by llndebetouw & Noriega-Crespol ( 11995b . According to this model, a bow 
shock seen on-axis at an inclination angle of ~30 degree, as in our case, shows qualitatively a 
behaviour like the one we observe, namely with the highest electron densities displaced at low 
velocity. 

In order to compute the contribution to the mass flux of the atomic component, we have de- 
rived M Kt from th e lum inosity of the [Fen] 1.644ytzm line, following the approach described in 



Garcia Lopez et al.l d2008l) . To compute the theoretical [Fen] lines fractional population, the n e 



values derived separately for the LV and HV components have been assumed. A constant value of 
10 000 K have been instead assumed for the electron temperature: the derived populations are only 
a little sensitive to this value and the results do not change significantly if the temperature is main- 
tained in between ~ 7000 K and 15 000 K. A n inclination angle of the jet with respect to the plane 
of the sky, i ~ 34° dEisloffel & Mun dt 1994), has been assumed and used to infer the velocity of the 
knots projected perpendicular to the line of sight through our measurements of the radial velocity. 
These values are reported in Table|6] while the adopted length of the knots projected perpendicular 
to the line of sight is listed in Table|5] The [Fen] 1.644 /xm line flux has been dereddened using the 
extinction derived from the ratio [Fe n] 1.644/1.257 and reported in Table[3] 

Mj et values have been derived both assuming a solar abundance of iron and considering the 
gas-phase Fe abundance measured from the Fe/P analysis, in the assumption that the same dust 
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HH34 HH 46-47 

Fell (1.64 Mm) Fell (1.64 |jm) 




-300 -200 -100 100 200 -200 200 



Velocity (km/s) Velocity (km/s) 

Fig. 8. PVD of the [Fe n] 1 .644 fan line for the HH34 and HH46-47 jets. 

depletion factor is valid for both the HVC and LVC. The corresponding range of values are listed 
in Table[6] Values of the order of 0.5-lxlCT 7 M yr _1 are found in the HVC, while the LVC shows 
~ 16 times lower values. 



7. Discussion 



The [Fe n] and H2 emiss ion structure at the base of HH46-47 is very similar to the one observed in 
several Class I jets (e.g.. 



Davis et al. 2001b 



2003 



Taka mi et al . 2006). As in the case of CTTS jets, 



the Forbidden-Emission Line (FEL) regions of Class I jets are characterised by the presence of two 
velocity components at the jet base at high and low velocity. The MHEL regions, as traced by H2 
emission, show on contrast just a LVC near the source while higher velocity H2 is usually detected 
further downsteam. In particular, the kinematics and physical parameters of the HH46-47 jet are 
very similar to those found in our previous work about HH34 dGarcia Lopez et al Loosl where the 
physical properties have been also derived as a function of velocity. For instance, Fig. [8]shows the 
[Fen] 1.644jt/m PVDs of the inner region of these two jets. Since the two objects are at a similar 
distance ( i.e. 450 pc), the emission regions have the same physical extent in this figure. In both the 
PVDs, the [Fen] lines broaden and emission at lower velocities, down to Okms -1 appears within 
~2" from the central source, i.e. up to distances of ~ 1000 AU. Th e extend of the LVC is comp arable 
with the typical size of the envelope of Class I sources (see, e.g., Andre & Montmerle1ll994 ). Thus 
this result may indicate an interaction of a low-velocity disc wind with the source envelope as the 
possible origin of the LVC in Class I jets. The large extend of the LVC is at variance with the PVDs 
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of most of th e CTT, where the LVC is usually confined within ~ 200 AU from the source dPvo et al 



2003 



2006) 



Also the values of electron density and mass flux rate that we derive in HH46-47, as well as 
their dependence with velocity are similar to those inferred on HH34, i.e., n e ~ 0.5-lxlO 4 cirT 3 in 



been also estimate d on a number of Class I YSO jets (e.g., 



Davis et al. 


2003; 


Podio et al. 


2006 



Takami et al. 2006). 



At variance with [Fen], the H2 PVDs of HH46-47 shows a different behaviour with respect 
to the one of HH34. In HH34 the H2 PVD shows spatially- and kinematically-separated LVC and 
HVC, and only the LVC is visible down to the central source. In HH46-47,on the contrary, a single 
velocity component near the source is detected. The morphology and kinematics of the H2 close 
to the source (the so-called Molecular Emission Line, MHE L, regions; Davis et alJboOlb ) in class 



I objects varies indeed significantly from source to source dDavis et al 



2001b 



2002) and this is 



another element that makes the origin of the MHEL regions unclear. 

In iGarcia Lopez et al.l (|2008) we have compared the [Fen] 1.644/im P VDs of HH34 with 



those of the synt hetic PVDs available for the disc-wind and X-wind models (IPesenti et al 



2004 



Shang et al.lll998l) . Most of this discussion is also valid for the HH46-47 jet. Briefly, none of these 



models are able to explain the persistency of the LVC at large distances from the launching regions. 
This may suggest that the low velocity that we observe far from the source represent dense gas en- 
trained by the high-velocity collimated jet, rather than gas directly ejected from the acceleration 
region. 

From the observed HH46-47 peak radial velocity of the HVC, corrected for the inclination 
angle, we estimate a jet poloidal velocity around 350kms _1 . This value is just at the border of the 
range of poloidal velocities predicted by X-wind models and it is also in agreement with disc-wind 
models with launch ing radii in the range 0.07-0.15 AU and magnetic lever arms larger than ~ 5 
( Ferreiraet ai]|2006 ). 

Further constraints on MHD launching models may be derived fro m the presence of dust in 

(2009). We have in fact seen that 



Podio et al 



the initial beam of the HH46-47 jet, as suggested by 
around 88% of the iron within 450 AU from the source is depleted, implying that a large fraction 
of dust grains are still present at the jet base. Dust grains of the ambient medium are expected 
to be destro yed by the p assage of the initial leading bow-shock, that have velocities higher than 
100 kms 1 ( Drainell2003 ). Since the amount of iron depletion decreases toward the external knots 
(only 58% of iron is found to be depleted in grains at ~5000AU from the source) dust grains 
are unlikely to come from entrained material. In addition they c annot be ref ormed along the jet 
beam since the time required for grain reformation is ~10 7 -10 8 yr ( JoneslboOl ). while the jet takes 
~1300yr to travel from the central source to the final bow-shock ( Hartigan et al.12005 ). Therefore, 
the larger amount of grains near the source may suggest that they have been removed from the 
accretion disc and transported through the jet. Interestingly it has been showed that, in the case of 
Herbig stars and TTauri stars, the disc is populated by dust gra ins only beyond the so-called dust 
evaporation point, R evp dAkeson et al. 2005; Eisner et al ]|2007blaL Future similar measurements of 



R evp for Class 0/1 sources would help to constraint the region from which dusty jets are launched 
from. 
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We have presented medium resolution (MR) (in the H- and K- band) and low resolution (LR) (in 
the J-, H- and K-band) spectra of the Class I driven HH46-47 jet. These observations allow us to 
resolve two velocity components at the jet base and to study the kinematics and physical properties 
of the jet as a function of velocity. From the [Fen] 1.644/1. 600 fim lines the electron density has 
been derived along the jet and as a function of velocity. In addition, from the luminosity of the 
[Fen] 1.644 /mi line the mass flux for the different velocity components has been derived. Finally, 
the [Fen] / [Pn] ratio has been used to retrieve important information about the dust depletion 
inside the jet beam. We summarise our results as follows: 

- We were able to trace the HH46-47 jet down to the central source. While the atomic component 
shows a HVC and LVC near the source, only a single velocity component has been detected for 
the molecular gas. Within 1" from the source, the atomic emission broaden down to Okms -1 
and even reaches red-shifted velocities. 

- From the ratio of the [Fen] 1.644, 1.257//m lines detected in the LR spectra, the visual ex- 
tinction for different knots have been derived. We have found an Ay ranging from ~6 mag to 
~12mag along the jet. 

- Several H2 lines have been detected from the LR spectra allowing us to construct excitation 
diagrams for different knots. From these diagrams we have derived the temperature and mass 
flux for the warm molecular component of the jet. The temperature has an average value around 
~2300 K, while we have found an average mass flux for the warm molecular component of 
> 3.7 x lO^Moyr" 1 . 

- The electron density, derived from the [Fen] 1.600/1. 644 /mi ratio, roughly decreases from the 
internal to the more external knots from 9400 to 3200 cm" 3 . Larger electron density values are 
found in the LVC with respect to the HVC for all the extracted knots, with an average value of 
-6000 and -4000 cm 3 for the LVC and HVC, respectively. 

- From the [Fen] / [Pn] ratio, we have estimated the iron depletion in the knots S and Zl. We 
have found that around 88% and 58% of the iron is depleted in grains in the knot S and Zl. 

- The mass flux rate transported by the [Fe 11] has been derived for both the velocity components 
assuming the Fe gas-phase abundance derived in this paper. We found that the HVC is trans- 
porting the bulk of the mass along the jet, with average valu es of M, et (HVC)~1.2xlO ~ 7 M yr~' 
and M ;e ,(LVC)~2xl0 _8 M o yr _1 . As in the case of HH34 ( Garcia Lopez et al^PPol . the differ- 



ence in tangential velocity between the HVC and LVC is not enough to account for the much 
hihger Mj et value in the HVC, indicating a significantly larger value of the (n H V |] product in 
the HVC than in the LVC (from a factor -1.5 in the knot Z2 to -3.5 in knot S). 
- Many of the derived properties of the HH46-47 are common to jets from YSOs in different 
evolutionary states. The derived densities and mass flux values are typical of Class I objects 
or very active T Tauri stars. However, the observed spatial extent of the LVC and the velocity 
dependence of the electron density are not explained by the current jet formation models. 
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2 where n H and V are the total density and the emitting volume, see iGarcia Lopez e t al. (2008) for more 
details. 
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Appendix A: knot fluxes at low resolution 



Table A.l. Observed unreddened lines along HH46-47: atomic lines. 



Line id. 






F (A F)(10 erg s 


1 -2\ 

cm ) 








c 


T TT T /I /r 1 DO 


LI 


A 


[Pn^Pz-'Dz 


1.188 


6.2(1.2) 








[Fen]a 4 D 5 / 2 -a 6 D 7/2 


1.249 


1.9(0.6) 








[Fen]a 4 D 7/2 -a 6 D 9/2 


1.257 


54.2(0.8) 


15.8(0.9) 11.9(0.7) 


7.1(0.9) 


6.4(1.0) 


[Fen] a 2 G9/2-a 4 D 7 / 2 


1.267 


3.2(0.8) 








[Fen]a 4 D 1/2 -a 6 D 1/2 


1.271 


2.0(0.3) 








[Fen]a 4 D 3/2 -a 6 D 3/2 


1.279 


2.4(0.6) 








[Fen]a 4 D5/2-a 6 D 5 /2 


1.295 


8.9(1.1) 








[Feii]a 4 D 3/ 2-a 6 D 1/2 


1.298 


3.4(0.9) 








[Fen] a 2 G7/2-a 4 D3/2 


1.300 


1.7(0.6) 








[Fen]a 4 D 7 /2-a 6 D 7 /2 


1.321 


13.1(0.4) 


4.5(1.4) 3.5(0.5) 


3.0(1.2) 


2.4(1.2) 


[Feii]a 4 D5/2-a 6 D 3 /2 


1.328 


4.3(0.6) 








[Fen] a 4 D 5 / 2 -a 4 F9/ 2 


1.534 


9.9(1.0) 




0.7(0.4) 




[Fen] a 4 D 3/2 -a 4 F 7/2 


1.600 


11.2(1.0) 




0.9(0.4) 




[Fen] a 4 D 7 / 2 -a 4 F9/ 2 


1.644 


85.8(0.7) 


8.9(0.5) 8.7(0.8) 


20.1(0.5) 


4.4(0.5) 


[Fen] a 4 Di /2 -a 4 F 5/2 


1.664 


5.7(0.8) 




1.1(0.7) 




[Fen] a 4 D 5/2 -a 4 F 7/2 


1.677 


13.3(0.9) 




1.8(0.4) 




[Fen] a 4 D 3/2 -a 4 F 3/2 


1.797 


5.2(0.9) 








[Fen] a 4 D 5/2 -a 4 F 5/2 


1.800 


14.3(0.7) 








[Fen] a 4 D 7/2 -a 4 F 7/2 


1.810 


52.6(1.1) 








[Fe n] a 2 P 3/ 2-a 4 P 3/2 


2.133 


3.2(0.5) 









Table A.2. Observed unreddened lines along HH46-47: H2 lines. 



Line id. 


A(/j.m) 






F(AF)(10- 


16 erg s 'cm" 


- 2 ) 








S 


Z2 


X 


Y3 


Y2 


Yl 


H 2 2-1 S(5) 


1.945 












2.7(0.7) 


H 2 1-0 S(2) 


2.034 




2.1(0.6) 


26.6(0.6) 


10.0(0.6) 


13.2(0.4) 


8.7(0.4) 


H 2 2-1 S(3) 


2.074 






3.7(0.5) 


2.3(0.2) 


3.9(0.4) 




H 2 1-0S(1) 


2.122 


12.0(0.5) 


5.6(0.3) 


36.3(0.3) 


27.0(0.3) 


3.3(0.2) 


22.7(0.2) 


H 2 2-l S(2) 


2.154 






1.1(0.2) 




1.6(0.2) 




H 2 1-0S(0)1 


2.224 


12.7(0.6) 


3.1(0.6) 


7.7(0.4) 


5.4(0.5) 


6.9(0.3) 


5.2(0.2) 


H 2 2-1 S(l) 


2.248 




1.4(0.4) 


3.4(0.3) 


3.0(0.5) 


4.6(0.2) 


3.5(0.3) 


H 2 1-0 Q(l) 


2.407 


16.6(1.3) 


4.8(0.6) 


28.9(0.6) 


20.8(0.8) 


24.3(0.7) 


17.8(0.5) 


H 2 1-0 Q(2) 


2.413 


13.8(1.5) 


2.3(0.3) 










H 2 1-0 Q(3) 


2.424 


16.2(1.3) 


5.8(0.5) 


30.7(0.6) 


23.3(1.5) 


26.5(0.7) 


17.8(0.5) 


H 2 1-0 Q(4) 


2.437 


9.3(1.1) 




8.4(0.7) 


6.1(1.2) 


8.0(0.6) 


5.2(0.6) 


H 2 1-0 Q(5) 


2.455 






23.4(0.9) 


15.7(1.6) 


21.6(0.8) 


14.9(0.7) 


H 2 1-0 Q(6) 


2.476 










6.6(0.6) 





11 Blend with [Fen] b 2 F 5/ 2-a 2 F 7/ 2 
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Appendix B: Rovibrational diagrams 



rot Y3 




6 8 10 12 14 

E v J (10 3 K) 

Fig. B.l. As Fig.|6]for knot Y3. 



rot Y2 




(10 3 K) 



Fig. B.2. As Fig.|6]for knot Y2 
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6 8 10 12 



E v J (lO^ K) 

Fig. B.3. As Fig.[6]for knot Yl 




6 8 10 12 



E v J (lO^ K) 

Fig. B.4. As Fig.|6]for knot X 
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